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Abstract 

In this work we derive local gradient and Laplacian estimates of the Aronson- 
Benilan and Li-Yau type for positive solutions of porous medium equations posed on 
Riemannian manifolds with a lower Ricci curvature bound. We also prove similar 
results for some fast diffusion equations. Inspired by Perelman's work we discover 
some new entropy formulae for these equations. 

1 Introduction 

The porous medium equation (PME for short) 
(1.1) atu = Au", 

where m > 1, is a nonUncar version of the classical heat equation (case m = I). For various 
values of m > 1 it has arisen in different applications to model diffusive phenomena like 
groundwater infiltration (Boussinesq's model, 1903, with m — 2), flow of gas in porous me- 
dia (Leibenzon-Muskat model, m > 2), heat radiation in plasmas (to > 4). liquid thin films 
moving under gravity (m = 4). crowd-avoiding population diffusion (to = 2), and others. 
The mathematical theory started in the 1950's and got momentum in recent decades as a 
nonlinear diffusion problem with interesting geometrical aspects (free boundaries) and pe- 
culiar functional analysis (like generating a contraction semigroup in and in Wasserstein 
metrics) . We refer to the monograph [V4j for an account of the rather complete theory con- 
cerning existence, uniqueness, regularity and asymptotic behavior of PME, mostly in the 
setting of the Euclidean space and on open subsets of it, as well as the different applications. 

The mathematical treatment of PME can be done in a more or less unified way for all 
parameters m > 1. Our main estimates below are only valid for nonnegative solutions, hence 
we will keep the restriction u > 0. This is reasonable from physical grounds since u represents 
a density, a concentration, a temperature or a height in the usual applications. However, re- 
writing (jl.ip in the more general form dtu = Aduj™"^?/), solutions with changing sign can 
also be considered, but the theory is less advanced. It has been proved that for given initial 
data uq S L^(]R") with uq > 0, there exists a unique continuous weak solution u{x,t) > 
of the initial value problem of (jl.ip . with a number of properties. 



Some of the existence, uniqueness and regularity properties hold true for the so-called 
fast diffusion equation (FDE), which is equation (jl.ip with m £ (0,1). FDE appears in 
plasma physics and in geometric flows such as the Ricci flow on surfaces and the Yamabe 
flow. However, there are marked differences between PME and FDE that justify a separate 
treatment of FDE, cf. |DK| . |V3j . In particular, the qualitative properties of FDE become 
increasingly complex for small m, far away from m ~ 1 (very fast diffusion). 

As is typical of nonlinear problems, the mathematical theory of PME and FDE is based 
on a priori estimates. In 1979, Aronson and Benilan obtained a celebrated second-order 
differential inequality of the form f |AB| ) 



(1-2 > TT- Sr- > , n ■■= 7^ 

^ ' ^ dx^\ dx' J - t n{m - 1) + 2 

which applies to all positive smooth solutions of (II. ip deflncd on the whole Euclidean spaccQ 
on the condition that m > 1 - 2/77,. Note that J2 ^ (mw™-^^) = A (^:^u"'~^ 

when TO 7^ 1. Precisely for the heat equation formula (|1.2p takes the form 

(1.3) Alogu+->0. 

obtained by setting m = 1. Estimate (|1.2p has turned out to be a key estimate in the 
development of the theory of the PME and FDE posed in the whole Euclidean space. In 
|AB] the estimate was used to prove the existence of initial value problem for PME and 
FDE. However, it has turned out difficult to find variants of (|1.2p that hold for flows posed 
on open domains, unless in dimension n = 1. 

In 1986 Li and Yau studied a heat type flow on manifolds ( [LYj ) . Among other things, 
they proved the following Li- Yau differential Harnack inequality. If {M",g) is a complete 
Riemannian manifold with nonnegative Ricci curvature and u:Mx[0,oo)— s-Risa positive 
solution to the heat equation dtu = Au, then there is a lower bound for Ag \ogu that has the 
precise sharp form (|1.3p . This extends to the manifold setting the Euclidean case mentioned 
above. Li and Yau also proved a local result, which implies (jl.3p when m is a global positive 
solution and M is a complete Riemannian manifold with nonnegative Ricci curvature. More 
precisely, they proved the following theorem. 

Theorem 1.1 Let {M",g) be a complete Riemannian manifold satisfying Ricci curvature 
Ric(Af) > —K^ for some K > 0. Let B{0, 2R) be a ball of radius 2R centered at O. Assume 
that u{x,t) is a positive smooth solution to the heat equation on B{0,2R) x [0,oo). Then 
for any a > 1 the following estimate holds on B{0, R) 

sup -a—]< ( — + KR^ ■ 

B(0,R) \ " 



uy - i?2 - 1 y 2(a- 1) 2t 



Here C is a constant only depending on n. 

The theory of PME and FDE on manifolds has not been considered until recently. The ex- 
tension of the Aronson-Benilan estimate (|1.2p to the PME on a closed Riemannian manifold 
with nonnegative Ricci curvature was done in the book |V41 Chapter 10]. 



^In dimension n = 1 the restriction is m > for general solutions, but we may keep m > — 1 for so-called 
maximal solutions ( [ERV] '). 
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In this paper we prove an extension of the Aronson-Bcnilan estimate to the PME flow 
for all m > 1 (Theorem 13. ip and the FDE flow for m E {rric, 1) (Theorem 14. ip on complete 
Riemannian manifolds with Ricci curvature bounded below. The estimates are of local type, 
hence even on Euclidean space, they give more information. The estimates look much better 
when the Ricci curvature is nonnegative. 

Recall that for the positive solution u ^ ^„ , such that € M^^'^(Af), to the heat 
equation, it was shown in [Nl] that 



(1.4) 



dW 
~dt' 



2t 



M 



where 



W(t)= / {t\Vf\^ + f-n)udti. 
Jm 



Using a basic identity involved in the proof of Theorem l3.11 we also obtain entropy formulae 
in the style of Perelman ([P]). This new entropy formula is the PME/FDE analogue of 

Organization. In §2 we introduce the main ideas of the regularity question for the PME 
in the Euclidean setting. We then pose the problems that have to be addressed. ^Sl contains 
the new local estimate for positive solutions of the PME on Riemannian manifolds with Ricci 
curvature bounded below. The estimate admits a version valid for the FDE if to € (toc, 1), 
which is developed in 21 Consequences in the form of Harnack inequalities for PME and 
FDE are derived in Sections [3] and [4] respectively. [jS] introduces and studies the entropies. 



2 Regularity of solutions of PME and Aronson-Benilan 
estimate 

A key idea in the PME theory comes from the observation that we can write the equation 
as a diffusion equation for a substance with density u{x, t) > 0: 

(2.1) Qtu = V • (c(u)Vw), 

We find a case of density-dependent diffusivity, i.e., c{u) = mu"^~^, so that c vanishes at 
u = 0; this makes the equation degenerate parabolic. It also implies the property of finite 
propagation, appearance of free boundaries, and limited regularity. Typical solutions with 
free boundaries are only Holder continuous in space and time. 

The second key idea in the study of the PME is to write the equation as a law of mass 
conservation 

(2.2) 5tu = -V-(uV), 

which identifies the speed as V = — tou'"^^Vw, and this in turn allows to write V as a 
potential flow, V = — Vp. This gives for the potential the expression p = mu"''-~^ / {m — 1). 
In the application to gases in porous media the potential is just the pressure and the linear 
speed-pressure relation is known as Darcy's law. Historically, the letter v has been used 
for the pressure instead of p, and we will keep that tradition. This variable v is crucial in 
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the study of free boundaries and regularity. Note that the pressure v := mu^'^ ^ /{i^ ~ 1) 
satisfies 

(2.3) dtv ^ {m-l)v/\v + \Vv\^. 

Thus, near the level u = we have the formal approximation dtv ^ |Vup, that can be 
easily identified as movement of the front with speed —Vv. It also means that the equation 
is approximately first-order so that we expect the Lipschitz continuity of v near the free 
boundary. 

Now we turn to regularity estimates for PME. The question of Lipschitz regularity of the 
pressure was solved in one space dimension, n = 1 , by Aronson who proved a local estimate 
for Vx using the Bernstein technique ( | Ar j l : a bounded solution defined in a cylinder in 
space-time Q = [a,b] x [0, T] has a uniform bound for \vx\ inside the domain, i.e., in 
Q' = [a', 6'] X [T',r], with a < a' < 6' < 6, < T' < T. Bf§nilan proved that in that 
situation vt is locally bounded in a similar way ( |Be| ). 

The extension of such results to dimension n> 1 fails, even for globally defined solutions. 
Indeed, it was shown in 1993 that the so-called focusing solutions are not Lipschitz continu- 
ous at the focusing point ([AG]), though wide classes of solutions can be Lipschitz continu- 
ous under special conditions ( |C VW| ) . And other kind of pointwise gradient estimates also 
failed. The problem of minimal regularity was reduced to proving Holder regularity, and 
this was done around 1980 by Caffarelli and Friedman ( |CF1| . [CF2j ). The proof of these 
results and the whole theory of the porous medium equation in several space dimensions 
was greatly affected by the existence of special one-sided estimates that we discuss next. 

The Aronson-Benilan estimate (11.21) can be written as 



(2.4) Aw > — — , when m > TOc, 

where we define v = logu for m = 1. Note that with this definition v <Q for < ?ti < 1 so 
care must be taken in manipulating inequalities when dealing with fast diffusion. Using the 
pressure equation (|2.3p . it immediately implies that for PME with to > 1 

(2.5) Vt > |Vvp - ^ ^v. 

so in particular 

(2.6) Vt > —(to — l)Kv/t, and ut > —Ku/t. 

Other forms of parabolic Harnack inequalities follow from such estimates, and lead to Holder 
regularity statements easily. These estimates have been used for all kinds of purposes in 
the theory, like existence of solutions in optimal classes of data, or asymptotic behavior, cf. 
IVl. 



A striking property of the Aronson-Benilan estimate is the fact that the constant k is 
optimal when m > TOc. Recall the Barenblatt solutions, in terms of the pressure, is given 

(2.7) = (£i!:^l_lff!)±, OO. 
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Equality holds in (|2.4p for v = Vc on the set {Vc > 0}. When m = 1 the estimates (|1.2p is 
optimal since equality holds in (|1.3[) for the Gaussian kernel. In some sense the Barenblatt 
solutions play for the PME a role that the Gaussian kernel plays for the heat equation. 

Many attempts have been made to obtain an extension of the estimate or a suitable variant 
for problems where the PME or the FDE are not posed on the whole Euclidean space: this 
can take the form of boundary value problems in bounded domains of M.'^, the PME posed 
on a Riemannian manifold, or even better, a local estimates valid in one of the two above 
settings. A straight extension of the global estimate to boundary value problems in several 
dimensions has not been done. In the case of homogeneous Dirichlct problems a literal 
extension is even false, in view of explicit solutions. So there is a hope for local estimates. 
In one space dimension, the local estimate of Vxx, hence of Vt, from below was obtained 
in |V1| , and the bound has a correction term involving distance to the boundary. But the 
method fails for n > 1 because it uses the previous knowledge of the local bound oi v^- A 
modified version of the local estimate will be the first objective of the present paper. 

Another research direction concerns the extension of Aronson-Benilan estimate to other 
equations, like the p-Laplacian equation or reaction-diffusion equations. Some work have 
been done, for example, in |EV| for the p-Laplacian heat equation on Euclidean space, and 
recently in [KNj for the p-Lapacian heat equation and (local and global) doubly nonlinear 
equation on manifolds. 



3 Local Aronson-Benilan estimates for the PME 

We proceed now with the new estimates. Let u > be a solution to the Porous Medium 
Equation (jl.ip . m > 1, posed on an n-dimensional complete Riemannian manifold (M",^). 
We will assume at least a local bound from below for the Ricci tensor. The initial and 
boundary-value problems for this equation are usually formulated in terms of weak solutions, 
or better continuous weak solutions |V4| . Our local estimate is closer related to the result 
of Li-Yau mentioned in the introduction, than that of say [SZ| . 

We will work with the pressure v, which satisfies equation (|2.3p . We see that Vv = 
mu'^~^'Vu, and in the case m < 2 this equation only makes sense over u > 0. In order 
to avoid this and other regularity difficulties in our computations we will assume that the 
solutions are positive and smooth everywhere. The smoothness property comes from lo- 
cal boundedness and positivity of u in view of standard non-degenerate parabolic theory. 
Application of our results for general weak solutions proceeds in a standard way by ap- 
proximating, and using the maximum principle and the local compactness of the classes of 
solutions involved. We refrain from more details on this issue, cf. |V4j . 

[Sjl Assuming that u > we introduce the quantities y = |Vup/u, and z = vt/v and the 
differential operator 

d 

C:= — -{m- l)vA. 

We also introduce the differential expression '■= az — y. Using equation (j2.3p we can 
write the equivalent formulae 

Vf iVwP 
(3.1) Fa = {m- l)Av +{a- 1)— = a(m - l)Av + {a - I)- -. 
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In particular, Fi = (m — l)Aw. Though our main goal is to estimate Fi, we will use the 
localization technique of Li and Yau to estimate F^ for a > 1. One reason is that sometimes 
the estimate of Fi is not feasible, for instance, when we want to obtain local estimates. 

The goal of this subsection is to calculate a formula for C{Fa). The following formula is 
helpful in the calculation: 

(3.2) c(^^^ ^-^C{f)-j^C{g)^2{m~l)v(v(^^^,\I\ogg 
Direct calculation shows the following Bochner-type formulae: 

Lemma 3.1 Let u be a positive smooth solution to on manifold {M"',g) for some 

m > 0, and let v :~ ^/jg pressurfl Then we have 

(3.3) Civt) = 2{Vv,\/vt) + Fivt, 
£(|Vwn = 2|Vz;pFi + 2(V(|Vwn,V'y) 

(3.4) -2(m - l)vvfj - 2{m - l)vR,jV,Vj. 

Here we use the notations Vi instead of v^i to denote partial derivatives, and we write Vij to 
denote the Hessian tensor H{v)ij of v, while vf^ denotes the standard norm of the Hessian. 
Rij is the Ricci tensor and R^jV^Vj ~ ffic(Vw, Vw). 

The following proposition is a generalization of the computation carried out in Proposition 
11.12 of [V4]. 



Proposition 3.1 Let u and v be as in Lemma \3Tn Then 

£{Fa) = 2(m - l)vfj + 2(m - l)i?yf + 2m(Vi^„, Vw) 

(3.5) +^a-l)(^Y + Fl 



Proof. Using p.2p and Lemma [5TT] we have that 

jr(\^) = -{2\Vv\'F, + 2(S/{\S/v\'),Vv)) 

\ V J V ^ 

I V'L' I' 

-2(?7i - l)vf, - 2(m - l)R,jViVj - ^—^ 



+2(to - l)w ( V ( ) , V log w 



C(-) = -{2{Vv,Vvt)+FiVt) 

_!^^+2(m-l).(vn,Vlog. 



Putting together gives that 



C{Fa) = 2{m-l)v(\'Fa,\'\ogv) +2{m-l)vtj+2{m-l)R^jV^Vj 



+a — ri — 2 Li — a 1 ^ 

V V V V V 



^Recall that when m = 1, we interpret v = logu. 
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Using 

(Vv^VivF^)) = v{VvyF,,) + F^\Vv\^, 
we can rewrite the last term in the above formula for C{Fa) as 

-{Vv,aVvt - V(|V«n) = 2{VvyF,,)+2Fj-^. 

V V 

Hence we get 

C(Fa) = 2mv{\7Fa,V log v) + 2{m - l)vfj + 2{m - l)R,jV,Vj 

+a — ri — 2 bi — a 1 h 2l<a • 

V V V V V V 

Note the last five terms simplifies as 



az 



(z-y)- 2y{z - y) - azy + y^ + 2{az - y)y = (a - + (2 - yf- 



This completes the proof of the proposition. □ 
When a = 1, (|3.5p becomes the following formula in jV4| : 



(3.6) CFi = 2{m - 1)<^- + 2(m - l)R^JV^VJ + 2to(VFi, Vw) + F/. 

From this, for positive smooth solution u to with m > 1 on a closed Riemannian 

manifold of dimension n with nonnegative Ricci curvature, the following estimate follows 
easily from maximum principle (sec Proposition 11.12 of [V4| ) : 

„ (m — 1)k n 

(3.7) Fi > n 



n(rn — 1) + 2 

[3l2 Now we prove a new local estimate for PME on complete manifolds. We employ the 
localization technique of Li and Yau ( |LYj . see also [KNj ) . Denote by B{0,R) the ball of 
radius R> and centered O in {M"-,g), and denote by r{x) the distance function from O 
to X. 

Theorem 3.1 Let u be a positive smooth solution to PME il.l]) . m > 1, on the cylinder 
Q B{0,R) X [0,T]. Let v be the pressure and let w^^x '■— ^'^^^b{o.r)x[o.t] v. 

(1) Assume that Ricci curvature Ric > on B{0, R). Then, for any a > 1 we have 

(3.8) a—<aa 7 + ^^^(61+62(0:)) 

V V \t R'^ 

on Q' :~ B{0, R/2) x [0,T]. Here, a := „"y|^™Y)+2 ^ (™ ^ ^)i^7 '^^d, the positive constants 
Ci and 62(a) depend also on m and n. 

(2) Assume that Ric > — (n — 1)K^ on B{0, R) for some K > 0. Then, for any a > 1, we 
have that on Q' , 

(3.9) -a-< ac? (\ + Cz{a)K'v^^^ + aa^^^i [C^) + C{{I<R)) . 
V V \t J R'^ 
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Here, a and C2(a) are as before and the positive constants 6*3(0;) and C'[{KR) depend also 
on m and n. 

Acceptable values of the constants are: 

Ci :=40(to- l)(n + 2), ^2(0):= 



a - 1 

C3(a) ^ '\ C[{KR) := 40(m - 1)[3 + (n - 1)(1 + KR)]. 

a — 1 



Note that C[{0) = C_ 



Proof, (i) We start with an auxiliary calculation about suitable cutoff functions. We take 
a cut-off function rj{x) of the form rj{x) ;= 9 {r{x)/ R), where 9{t) is a smooth monotone 
function satisfying the following conditions 9{t) = 1 for < t < ^, 9{t) = for t > 1, 
{9'Y /9 < 40, and 9" > -409 > -40. 40 is just a convenient number, it could be optimized. 
Direct calculation shows that on B{0, R) 

and also 

f■,^U A ^ 40((n-l)(l + j^i?) + l) 2 

(3.11) At] > —5 , if Ric > -fn - IjK , 

R^ 

with the help of the Laplacian comparison theorem. In particular Arj > — 40ni?~^ when 
Ric > 0. We give some details about deriving (|3.1ip . Note that 



A77 : 



9" iVrl^ 9' At 



i?2 R ■ 

By the Laplacian comparison theorem when Ric > — (n — 1)K^ 

Ar <{n- 1)K coth{Kr). 
Since cotli is decreasing, and 9' = when r{x) < ^R, this implies 



-40 V40(n-1)^. wU.^N 

where we have used \9'\ < V40. Using the inequality K coth{KR) < -^(1 + KR), one has 

-40 2x/40(n-l), 



(ii) To obtain the desired estimates, we apply the operator C to the function tri(—Fa), 
then apply the maximum principle argument. Note that if tr]{—Fa) < on Q, then p.9p 
follows. Now we assume in.ax(^^ t^^Qtr]{—Fa) > 0. Let {xo,to) be a point where trj{—Fa) 
achieves the positive maximum. Clearly we have to > and at (xo,to) 

VF„ = CM-F^)) > 0. 
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All further calculations in this proof will be at {xq, to). Denote C4 := 40((n— + KR) + 
1), y := r/y = 77 ^^"^ and z := rjz = 77—. Combining (j3.5p with the above estimates of i], we 
have that when Ric > -(n - l)^^ on B{0, R), 

< T^CM-F^)) 



< - 



trf (2(m - l)v]j + 2(m - l)R,jV^Vj) + 2mtff {\I {-F^), Vw) - (a - l)trfz'^ 



^tifFl + 2t(m - \)v'q'\^^{-F^) + ■§(m - l)tvr^{-F^) + 7?2(_F«) 

TXiifHi — 1) I 

< ^ ^— ^ • t{y - z)^ + 2{n - l)(m - l)K'^ty-qv + 2mt{y - az)|Vr7||V7;| 

n(m — 1) 

-(a - + (y _ «5) (^(m - l)^^^^^^ -tv + iy 

In the last inequality above we have used v^j > ^^^^ and (m — l)Az; ^ z — y ^ Fi. Now 
write C5 := 80 + C4 and 

2 1/ n2 a — \ f ~ 2 

[y-z) =—{y^az) +2 — —[y-az)y+l y. 

\ a J 

Also note that 

2mt{y - q;z)|V?7| • |V-u| < —mt{y - az)y'^/'^v^/'^ . 

R 

Putting these together we have that when Ric > — (n — 1)K^ on i?(0, i?), 



< ttOj — az) +t{y — az 

^ ^ 2 

a \ a 



/ 2(a- 1). 40m ^^.^ . , , C5 \ 



(3.12) +{y-az)--(- -] + 2{m - l){n - l)K'^tyvri - {a - l)tz^ 



(1) When K = 0, using -Ax^ + < it follows from (pl^ 

tw / 200aa^rr 
'R^ \ a-l 



iA ' 

^^2 



t , ,0 , ^ / /200aa^m^ , , 
(3.13) < - — (y - + (y - a5) ( ^ ( — — + (m - 1)^5 ) + 1 



aa 

This gives the first estimate 



(2) When K ^0, in ([3^2]) we handle the {y - a5)-term as in ([333]) with Ce := ^""""i"" 
(m — 1)C5, and use 



1 / 1 \ 2 
1 / a — 1 



i?/^ + 2{m - l)(n - l)tK'^yvij < Cjtv'^ , 

■~ (a -1)2 ■ 

Then the above quadratical inequality (j3.12|) on (y — az) reduces to 

< *—{y - a~zf + f + 1 ) (^ - az) + Ct^'. 



a \ a 
where 
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aa^ / w 1 // V C7 



This implies that 



\ R'^ t a — 1 

The claimed result follows from this easily. □ 



Using the local estimate one can generalize Proposition 11.12 of |V4| to noncompact com- 
plete Riemannian manifolds with nonnegative Ricci curvature. 

Corollary 3.1 Let u{x,t), t e [0,T], be a smooth positive solution of the PME il.l]) with 
m > 1 on a complete manifold {M",g). 

(1) If{M,g) has nonnegative Ricci curvature, then Jg. 7\ l holds fort G (0,T], provided that 
v{x,t) = o(r^(a;)) uniformly in t £ (0,T]. 

(2) If the Ricci curvature Ric > — (71— 1)A'^ on M for some K >Q and Wmax ma-XMxfo.T] ''^ < 
00, then for any a > 1 

1^,,. vt |Vz;P ^fl , (m-l)(n-l) ^^, 

(3.14) a > -(to - l)Ka - H K w,„ax 

V V \t a ~ 1 

Proof. (1) Taking R 00 and then a ^ 1 in ()3.8p we have the result. 
(2) Taking i? ^ 00 in (|3.9[) we have the result. □ 

Integrating along minimizing geodesic paths of the local estimate, one can obtain the 
following Harnack inequality. Here we just state the most general form. When if = 0, if we 
assume Wmin minjv/x[o,T] v > 0, the estimate simplifies by taking a ^ 1. 

Corollary 3.2 Same notations and assumptions as in Theorem \3.1\ Denote v^J^'"'" to be 
miug^Q ii)x[o T] -Assume that Ric > — (n — 1)K^ on B{0, R) for some K > 0. Then for 
any xi,X2 G B{0, ^) and < ti < t2 < T , and any a > 1 

v{xi,ti) 

where d{xi,X2) is the distance and the constants C2(a), 6*3(0;) and C'i{KP) are as in The- 
oremWTTX 



Proof. For the minimizing geodesic ^[t) joining [xi,ti) and {x2,t2) we have 

,„,(4b44) . rvzi^fZi;,,)),, 



v{xi,ti)J Jti \v V 



av 4w 
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The result follows from the observation that 7(3) lies completely inside B{0, -j) and the 
estimate in Theorem 13. II □ 



A different way manipulating the integration on geodesic path can have the following 
consequence of Corollarv 13.11 This estimates are the analogue of the classical one for the 
positive solutions to the heat equation (in viewing for the heat equation v = log u) . 



Corollary 3.3 Same notations and assumptions as in Corollarv \3.1[ We further assume 
that Umax '■— iii^Xj\/x[o,T] ''^ < c«. Let xi, X2 S M and < ti < t2 < T . 

(1) If{M,g) has nonnegative Ricci curvature, then 

/ , N / + ^ ^ / 1^ 1 *2 d'^{xi,X2) 

v{X2,t2j -v{Xi,ti) > -(?71- IjKWniaxlog- -77 —- ■ 

ti 4:[t2 — ti) 

(2) // Ricci curvature Ric > —{n — 1)K'^ for some K > 0, then for any a > 1 

v{x2,t2) - v{xi,ti) 

t2 (m- 1)^(71- l)Ka ^ 2 /. , N a(P{xi,X2) 

> - (to - IjKaUmax log K Umax(^2 - h) - — — . 

ii a - 1 4(t2 - 1\) 

Proof. We only prove (2). Let ^{t) to be a constant speed geodesic with j(ti) = xi and 
7(^2) = X2- We compute using (|3.14[) 

/■t2 

v(x2,t2) - v{xi,ti) > / Vt + {'Vv,'j) dt 

Jtl 

> -(to - l)KaVmaxl0g- K V^^^[t2 - ti) - — — . 

ti a — I 4(t2 — tij 

□ 

Theorem 13. II can be used to give a local lower estimate for the Laplacian of fox (3 > 1. 

Corollary 3.4 Same assumption and same notations as in Theorem \8.1\ Let 1 < /3 < 
to/ {m — 1) fixed. Define a by = (to — l)(/3 — 1). 

(1) Assume Ric > on B{0, R), then we have on Q' 

(2) Assume that Ric > — (n — 1)K^ for some K > on B(0, R). Then we have on Q' 

Av^ > -nap {v^gf-' Q + C,{a)K'v^^£j 

-^aP^Jl^ (C2(a) + C[{KR)) . 
where the constants C2{a), Cs^a) and C((/\P) are as in Theorem lS.li 
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Proof. Clearly a > 1. We compute 



^ ^-(^-^( a{m-l)Av + {a~l)^— 



a{m — 1) 



a{m — 1) 

The corollary follows easily from Theorem 13. II □ 



4 Local Aronson-Benilan estimates for FDE 

The fast diffusfon equation, FDE, is equation (jl.ip with m e (0, 1). However, as we have 
seen, Aronson-Benilan estimate for FDE on Euclidean space holds only in the range 1 > 
m > rric '■= ^~ where the relevant constant k = n/ {n{m—l) + 2) is still a positive number. 
This is also the range where the Barenblatt solutions can be written and play similar role as 
they play in the theory of PME. Hence the Aronson-Benilan estimate holds on the range as 
it would be expected. There is another point needs to be made. Since m < 1, the pressure 
V = -^^^u™~^ is negative and moreover it is an inverse power of u. But (m — l)v is still 
positive, hence, as shown, the inequalities (|2.5p and (|2.6|) hold for all m > rric. 

[4ll Let u be a smooth positive solution to FDE p.ip with rn G {rric, 1) on a closed Rieman- 
nian manifold (M",g) with nonnegative Ricci curvature. From (|3.6p . and since m — 1 < 
we have 



£(^^1) - 2m(VFi, Vi>) <F^^ 2(1 - m)«f, < - ( " . - 1 ) F^ 

■' V (1 ^ m)n ' 

The following follows easily from maximum principle 

„ (ni — 1)k . K 

(4.1) ^^i<-^ ie., Az;>--. 



Note the direction in the Fi -inequality differs from ()3.7p because of m < 1; on the other 
side, the Aw-inequality is the same as in the case m > I. 

[4l2 Now we prove a new local estimate for FDE with m £ {nic, 1) on complete manifolds. 
In this subsection we employ the same notations as in §3.2 and use a similar localization 
technique as that of Li and Yau. It turns out that this case technically is slightly harder 
than the previous case. For example we have to make use of the term (a — l)tz'^ (which 
was simply dropped before). As we have seen from §4.1 for m £ {rric, 1) we should estimate 
Fa from the above (instead of from the below). For the local estimate, another difference is 
that for FDE we will estimate Fa for a < 1 (instead of a > 1). 

Let (xojio) be a point where function trjFa achieves the positive maximum. Clearly we 
have io > and at (a;o,to) 

\7Fa = -^Fa, C{tr,Fa)>0. 
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All further calculation in this proof will be at {xq, to). Combining (|3.5p with the estimates 
of 77, we have that when Ric > — (n — 1)K^ on B{0, R), 

< TjCitTjFo,) 

< trf (2(m - l)vfj + 2(m - l)R,jV,Vj) + 2mtr)'^(VFa,\'v) + {a - l)tr]'^z^ 

< -^ —^tiv - 5)^ - 2(n- 1)(7«- l)K'^tyriv + 2mt(az - y)|Vr;||Vt;| 

?t.(to — 1) 

+ M - y) (^(m - l)i^!^±^ . to + 1^ - (1 - a)tS^. 

Noticing that now we have (m — l)v > and a < for to e (yric, !)• Proceeding as in the 
proof of ((3?T2)) we then have that when Ric >-(?!- l)K^ on B{0,R), 



< - ^(az - y)^ + t{a~z - y) ^^^(-^) + ^ • + (m - 1)^ 

— aa^ \ — aa^ it rt^ 

+ {az - y) - — (^—^^ - 2{m - l)(?i - l)K^ty7]v - (1 - a)tz^. 



Here we can not do what was done in the proof of (|3.12p since ^i"''^^ "'^ the coefficient in 
front of —y (> 0) is positive in stead of being negative. For simplicity let X = az — y > 0, 
Y = —y >0, /?=1 — a>0, 7 = —a > 0. By writing = -^{ctz — y + y)^, and collecting 
terms we then have 

< ^— ({l + jf3)X^ ~2P{l+j)XY + p{l3 + -f)Y^) +2{n-l)[{m-l)v]K'^tY 

+ '-^XYH^vf^ + '§X[im-l)v]+X. 
Let v^'l^ := ma.XB{o.R)x[o.T]{~v)- Using that for any £1,62 > 

and 

we get 

< - ^ ((1 + 7/3)^' - 2/3(1 + 7 + ei)Xy + /?(/? + 7 - /^e^)^^) + x 
70;^ 

t7a^l600TO2 tCs^., 1)^(1 -TO)^(z;g,^)^g\ 

Choosing £2 small such that 
(4.2) P + ^-l3el>0, 
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then 

-2/3(1 + 7 + ei)XY + /?(/? + 7 - i3el)Y^ > -pSlJ^l^^x^ 
hence we have 



, (l+7 + ei)' 

(3 + j-l3el 



-/a^ \ " ' 13 + "f- Pel J 2ei/3 i?2 



Now for any (3 E (0, 1) we can choose small ei > and £2 > such that both ()4.2p and 

(4.3) l + 7/3-/3^^±^>0 

/3 + 7 - /3e^ 



hold. For example, any (ei, £2) S ^0, g^"l[° ^-^ x ^0, \J ^^^^ works. Repeating the argument 
m the proof of Theorem 13.11 after p.l2p . we have proved the following result for FDE. 

Theorem 4.1 Let B{0,R) he a hall in a complete Riemannian manifold (M"',g), and 
let u he a positive smooth solution to FDE 11. 1]) with m £ (tjc, 1) on the cylinder Q = 
B{0,R) X [0,T]. Let v he the pressure and let w^^x = b {o, R) x [0. t]{~ ">-')■ 

(1) Assume that Ric > on B{0, R). Then for any < a < 1, ei > satisfying 

(4.4) (7i(a, a, ei) 1 + (-a)(l - a) - (1 - ^±fil! > 0, 

(1 — a) — a 

we have 

(4-5 a—>—^ 7 7 + (<-2 a,a,ei) + C3) 

w w 6i(a, a, ei) \t R"^ ^ 

on Q' = B{0,R/2) x [0,T] w/iere 

C'2(a,a,ei) := 1600m^ „ , > 0, 
2ei(l - a) 

C3 :=40(l-TO)(n + 2) > 0. 

(2) Assume that Ric > — (n — 1)K^ on B{0, R) for some L{ > 0. Then for any < a < 1 
and £1 > 0, £2 > satisfying 1^4 .2^ and we have 



— "~ 5^ — \ t '-^4(a,£2)yCiK 

(4.6) - ''^f (C2(a,a,£i) + C.iKR)) 

on Q' . Here 

(l-a + £i)2 



C[ := C[{a, a, £1, £2) := 1 + (-a)(l - a) - (1 - a) 
C'4(a, £2) 



(l-a)-a-(l-a)£2 

(n — 1)(1 — m) 



(1 - a)£2 

C'5(A'i?) 40(1 - m)[3 + (71 - 1)(1 + KR)]. 
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Noticing that Ci(a, a, ei) ^ 1 as a — *■ 1. Similar to the proof of Corollary 13.11 by taking 
R ^ oo and then a — > 1, Theorem 14. II has the following consequence for global solution of 
FDE. 

Corollary 4.1 Let u[x,t), t e [0,r], be a smooth positive solution of the FDE with 
m £ {iTLc, 1) 071 a complete manifold {M",g). 

(1) If{M,g) has nonnegative Ricci curvature, then 

|Vwp Wt ^ (1 — m)K 



V V ~ t 

holds for t G (0,T], provided that \v\{x,t) = o(r^(x)) uniformly in t £ (0,r]. 

(2) If Ricci curvatureHic > —{n—l)K'^ on M for some K > andwmax t^^'^MxIo.t]{^'^') < 
oo, then for any a G (0, 1) and 62 > satisfying 

C'{ := C'{{a, a, £3) := 1 + (-a)(l - a) - (1 - a)- . " > 0, 

(1 — a) — a — (1 — ajej 



we have 



Vt^P vt ^Ko" (I , (l-m)(n-l) [^,^^.2- 
— a—> -(1 - m)-^ - + — — \^C'{K V, 



max 



V ■ ' C'{ \t (l-a)e2 

Integrating (|4.6p along minimal geodesic, we obtain 

Corollary 4.2 Same notations and assumptions as in Theorem \4-l\ Denote v^i^'"^ to be 
ming(-(»j n)x[o T]i^'^)- j4,ssume that Ric > — (n — 1)K^ on B{0, R) for some K > 0. Then 

6 ■ 



for any xi,X2 € B{0, ■^) and < ti < t2 < T , and any a > 1 



-v{x2,t2) ^ f t2\ ( ad'^{xi,X2) 

< — • exp ' 



-ViXl,tl) \tj \^yW^^^_^^^l 

■ exp [^^-^v^L [Ci{a,e2)^JC[K^ + ^ -) (^2 - ^0 

where the constants C'l, C2{a,a,ei), (74(0;, £2) and C^{KP) are defined as in Theorem \4.1\ 
Proof. For the minimizing geodesic j(t) joining {xi,ti) and (0:2, t2) we have 

\~v{xi,ti)J Jti \v V J 

The result follows from the observation that 7(3) lies completely inside B{0, -j) and the 
estimate in Theorem 14. II □ 



An integral version of Corollary 14.11 is 
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Corollary 4.3 Same notations and assumptions as in Corollary \4-l\ We further assume 
that Wmax < oo. Let Xi,X2 S M and < ii < t2 < T. 

(1) If{M,g) has nonnegative Ricci curvature, then 

/ , N / + \ ^ n \ - 1 *2 d'^{xi,X2) 

v{X2,t2) - v{Xi,ti) > -(1 - TOjKWniaxlog- -77 ■ 

ti 4:[t2 — ti) 

(2) // Ricci curvature Ric > — (n — 1)K^ for some K > 0, then for any a G (0, 1) 

V{x2,t2) - V{xi,ti) 

Ka _ , t2 (to — l)^(n — l)Ka 2-2 ad'^{xi,X2) 



> - (1 - TO)-^Umaxl0g- — = K W,„ax(i2 - ti) 



Corollary 4.4 Same assumption and same notations as in Theorem \4-l\ Let (3 G (l,oo) 
fixed. Define a e (0, 1) by ^ ^ (to - - 1). 

(1) Assume Ric > on B{0,R), then we have on Q' 

Li{a,a,ei) ^ ' \t R'' ^ 

(2) Assume that Ric > — (n — 1)K^ on B{0, R) for some K > 0. Then we have on Q' 



^1 V - max/ ^ ^ 

+ !^(!^(C2(a,a,ei) + (75(i^i?)). 

5 Entropy formulae 

In this section we show that the computation in fJ3] and thermodynamic considerations from 
Perelman's [P] (see also |N11 IN2| ') lead to a formula like Perelman entropy formula for the 
porous medium equation. The same formula works for the fast diffusion equation and with 
good concavity for 1 — i < to < 1. 

Assume that M is a compact manifold. First we derive some auxiliary integral formulae. 

Lemma 5.1 Let u be a positive smooth solution of with to > 0, and let v, Fi be as in 

gl Then 

(5.1) — / vudfi = I Fivud^^—m f jVupud/z. 

dt Jm Jm J m 

Proof. By (|l.ip and (|2.3[) wc have dt{vu) ~ {m ~ l)vuAv + |Vwpu + v/S.u^ . Recall now 
that Fi = {m - l)Av. Then, 



d 
di 



/ vud^i= / Fivud^i+ / (IVwpw + wAu'") d/i. 

JM JM JM 
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Using the identity Vm™ = uVw, a simple integration by parts shows that the second term 
on the right-hand-side is zero. The first part (|5.ip foUows. Integration by parts shows that 



□ 



/ Fivudp= I {m — l){Av)vudfi = m / (Au)m'" d/i = — m / iVwpwd/i. 

JM JM JM JM 



Lemma 5.2 

(5.2) ^ / Fivudfi = 2/ {{m- I) {vfj + R^jViVj) + F^) vudfj,. 

J M JM 



Proof. Using (jl.ip . (|2.3p . and the formula dt ^ C + [m — l)wA, we have 

d 
di 



Fivudfi = / dtFivu + Fidt{vu) dfi 

M JM 



{CFi)vudii+ / {m — l)v{AFi) vud^ 

M JM 

Fi[{m - VjvuAv + |Vwpw + uAu™] d^i. 

IM 

We will use p.6p to compute the term with £, and we also use (m — l)Aw = Fi. Then, 

/" Fivud^i = [ {{m- l)v{AFi)vu + F^vu + FivAu"") d^i 
dt Jm Jm 

-|-2(to— 1) / {vfj+RijViVj)vudfi+ / Ffvudji 

JM JM 



+2m / {VFi,Vv)vudn+ / Fl\Vv\'^ud^l. 

JM JM 



IM JM 

Using identity (m — I)V(w^m) = {2m — l)vu\/v, and integrating by parts, we get 



(m-1) / {AFi)v^ud^i^ -{2m-l) / {VFi,Vv)vudfj,. 

JM JM 

Finally, since Vw™ ~ uWv we also have 

/ FivAvJ^'dpL = - / (VFi,Vu)vudAi- / Fi|Vwpwd/i. 

JM JM JM 

Combining these equalities we prove (|5.2[) . □ 
Recall the constant a — "^™7,^? „ = (m — 1)k. Wc first define 

n(?n— 1)+2 ^ 

(5.3) Nu{t) ~ -t" I vud^i. 

JM 

By Lemma |5.1[ we have that 

d 



(5.4) -Mu{t) = (^F, + ^) vu dfi 
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Note that the universal estimate (|3.7p amounts to Fi + | > 0. Now we define the entropy 
to be 



(5.5) 



Using the last part of ()5.ip to compute the term with Fi , we get 

|Vw|2 a+1 



Wu{t) = t" 



m- 

M \ V 



t 



vu dji. 



Now we show the following Perelman type entropy formula for PME. We put b = n{m— 1), 



so that a — 



Theorem 5.1 Let u be a positive smooth solution to lil.l]) with m > 0. Let v be the pressure 
and let Wu(t) be the entropy defined above. Then 



di 



(5.6) 



a+l 



M 



{b + 2)t 



9i3 



+ RijViVj \ vudu 



-2r+i / (Fi + -) vudn. 



Proof. Note 



By Lemma 15.11 it is easy to see that 

Jm 



dt V dt 



Hence, 



dt 



Wuit) 



dt 



M 



((m - 1) + RijViVj) + F1) vudfi^ 



M 







1 (FH 




J M ^ 





-2t 



a+l 



M 



((to - 1) [v^j + RijV^Vj) + FI) vudfi 



-(a + l)i"/ (Fi + -)vudfi- {a + l)^ Fivudyi 
JM ^ Jm 

(to - l)vf, + (to - l)^(Ai))^ + (to - 1) '^^^ Av 
iM \ t 



-2t 



a+l 



a-^ + a 
~^t^ 



vu dfi — 2{m — l)t°' I RijViVjVudfi. 

M 
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Observing that a + 1 = ^ = = (m - 1)^^, hence 

(m - \)v% + (m - \f{l\vf + (m - l)^^Az; 



= (m - 1)4 + (m - l)^(At-)^ + 2(m - 1) (^:p^A« + (m - 1)- ^ ^ 2)2^2 
/ 2 2 n 



(6 + 2)i (6 + 2)2t2 

Formula (j5.6|) follows from completing the squares. □ 



Corollary 5.1 Let {M,g) be a closed Riemannian manifold with nonnegative Ricci curva- 
ture. Assume that u is a positive smooth solution to PME il.l]) with m > 1. Then 

(1) ^Muit) < and jf.yVu{t) < 0. In particular Afu{t) is a monotone non- decreasing 
concave function in i. 

(2) Any ancient positive solution to must be a constant. 

Proof. We only need to justify (2). For the ancient solution, (|3.7p implies Fi > 0. On the 
other hand, 

/ Fiuvdfj, — —m / \Vv\'^udfi 
Jm jm 

which implies that |Vw| =0. □ 

Remark 5.2 The result can be proved for complete noncompact Riemannian manifolds with 
the help of the gradient estimates from the previous section. The interested readers can find 
the details for the linear heat equation case in \Chf . 



In |CTj . the following logarithmic Sobolev type inequality related to PME was proved on 
S". (See also [DD| for Sobolev type inequalities related to FDE.) 

For any m > 1, / G ^^(R") n i™(M"), |V/'"-5| e L^{W"), we have 

2 



1 \ r < y 2m 



m — 1 / /»„ 2 V 2m — 1 



vr 



dx + A,n{\\f\\l), 



where Am{K) ■= /r„ (^^-^oo + :;^u^^ dx with Uoc ■= (C - ^l;;rl^l^)+ ' ^^'^'^9 i^e 
Barenhlatt-P attic solution of order m and mass K := J^,^ Uoo dx. 

The previous defined entropy yV„ is related to the above Sobolev inequality in the following 
way. First notice that A(m, n) := A„i{l) is a constant depending only on the dimension 
n and m. Direct calculation shows that Carrillo-Toscani's logarithmic Sobolev inequality 
amounts to that 

I \ . / 1 



For the fast diffusion equation, the monotonicity of entropy W„ is similar to that of PME 
given in Corollary 15.11 
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Corollary 5.2 Let [M,g) he a closed Riemannian manifold with nonnegative Ricci curva- 
ture. Assume that u is a positive smooth solution to FDE il.l]) with m < 1. Then 

(1) 4iNu{t)<Qforme (m„l). 

(2) -^iV^uit) < for m G [m'^,1) with m'^ = 1 — i. In particular Afu{t) is a monotone 
non- decreasing concave function in j when m G [to^, 1). 

(3) Any ancient positive solution to with m G {mc, 1) must be a constant. 



Proof. (1) In dnH), t; < and by Corollary SU^l) Fi + f < 0. 
(2) Notice that 

2 



-(to — 1) 



1 



ib + 2)t 



9i] 



> 



1 



n(m — 1) 
1 



(to, — l)AiJ + 



n{m — 1) 
{b + 2)t 



Fi 



n{m — 1) 

Note V < 0, the result then follows from (|5.6p and 1 + (to — l)n > for to G [toJ,, 1) 



□ 
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